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ABSTRACT 


During  May,  June  and  July  of  1 965  the  received  field 
intensity  frapn  both  VLF  stations  on  Oahu  Island,  Hawaii  were 
measured  at  Tokyo,  Japan;  Guam;  and  Boulder,  Colorado. 

Atmospheric  noise  field  intensities  were  also  measured  at  all 
three  sites  and  relative  phase  variations  of  the  received  signals 
were  measured  at  Japan  and  Guam.  This  program  presented  a 
unique  opportunity  to  study  VLF  propagation,  since  the  paths 
from  Oahu  to  Guam  and  Japan  were  almost  the  same  length  and 
were  both  all  sea  paths.  The  availability  of  two  transmitters  on 
Oahu  which  permitted  simultaneous  measurements  at  two  VLF 
frequencies  over  the  same  path  was  also  a  unique  opportunity. 

The  results  of  this  measurement  program  are  presented  in  a 
number  of  ways  in  this  report  to  afford  maximum  opportunity  to 
study  the  propagation  and  atmospheric  noise  conditions. 

The  results  show  22.  3  and  24.  0  kHs  to  be  the  best  frequencies 
of  those  measured  for  the  Oahu  to  Guam  and  Japan  paths.  Significantly, 
the  diurnal  variations  of  field  intensity  were  very  repeatable  even 
during  sunset  and  sunrise  transition  times  which  indicates  that 
better  predictions  of  VLF  field  intensities  may  be  possible  with 
an  improved  propagation  model. 

Means,  medians,  standard  deviations,  cumulative  probabilities 
and  diurnal  variations  of  signal,  noise  and  signal-to-noise-ratios 
have  been  computed,  plotted  and  tabulated  for  all  frequencies  and 
for  four  hour  time  blocks  and  twenty-four  hour  periods.  In  addition, 
continuous  plots  of  signal  data  for  Japan  and  Guam  are  presented. 


1.  INTRODUCTION 


Although  many  advances  have  recently  been  made  in  the  science  and 
art  of  predicting  VLF  propagation  characteristics  and  atmospheric  noise 
levels,  there  still  exist  unknown  factors  and  discrepancies  between  predictions 
and  results.  The  refinement  of  prediction  techniques  calls  for  refined  field 
measurements  to  attempt  to  resolve  some  of  the  problems.  This  report  presents 
the  results  of  such  a  measurement  program  during  May,  June  and  July  of  1965. 
The  program,  equipment,  field  operations,  final  results  and  analysis  of  results 
are  presented  herein. 

2.  DESCRIPTION  OF  PROGRAM  AND  EQUIPMENT 
2. 1  Transmitter  and  Receiver  Locations 

Figure  1  shows  a  map  of  the  World  indicating  the  location  of  transmitters 
and  receiver  sites.  The  location  of  two  transmitting  facilities  on  Oahu  Island 
afforded  the  unique  opportunity  to  measure  propagation  characteristics  along 
identical  paths  at  two  frequencies.  Receiving  sites  located  both  west  and  east 
of  the  transmitters  offered  opportunity  to  determine  directional  and  latitude 
effects. 

2.  2  Frequencies  and  Parameters  Measured 

Frequencies  of  16.  6,  19.  8,  22.  3,  24.  0  and  26. 1  kHz  were  used 
at  various  times  throughout  the  tests.  The  Haiku  transmitting  facility 
was  primarily  operated  at  19.  8  kHz  as  a  means  of  providing  a  reference 
transmission  for  better  delineation  of  effects  seen  at  other  frequencies. 

Some  transmissions  at  16.  6  kHz  were  made  from  Haiku  on  a  short  term 
basis.  The  Lualualei  transmitter  facility  transmitted  on  all  of  the  frequencies 
at  various  times  during  the  tests. 

The  receiving  sites  at  Guam  and  Japan  were  equipped  with  identical 
equipment  designed  to  measure  the  following: 
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1.  The  vertical  electric  field  intensity  from  both  transmitters. 

2.  Phase  variations  of  the  arriving  signal  from  both  transmitters. 

3.  The  average  vertical  electric  field  intensity  of  the  atmospheric 
noise  at  the  frequencies  of  transmission. 

4.  Semi-automatic  recording  of  atmospheric  noise  field  intensities 
at  frequencies  from  30  Hz  to  30  kHz. 

5.  Amplitude  probability  distributions  of  atmospheric  noise  at 
frequencies  from  50  Hz  to  30  kHz. 

6.  Direction  of  arrival  of  atmospheric  noise. 

At  Boulder  equipment  was  available  only  for  measuring  the  vertical 
electric  field  intensity  of  the  arriving  signals  and  of  the  atmospheric  noise 
at  the  frequency  of  transmission. 

2.  3  Equipment  Description 

Two  Tracor  VLF  receivers  were  shipped  to  both  Guam  and  Japan  for 
amplitude  and  phase  measurements.  Ten  foot  whips  and  associated  couplers 
were  used  with  these  receivers  for  measuring  the  vertical  electric  field. 
Rubidium  frequency  standards  were  used  for  the  reference  frequency  of  these 
instruments.  Hewlett  Packard  302-A  Wave  Analyzers  were  also  shipped  to 
each  station  as  a  back  up  instrument  for  obtaining  vertical  electric  field 
intensities  and  for  atmospheric  noise  measurements.  These  instruments 
were  to  be  used  with  30  foot  whips  and 'their  associated  coupler.  Instrumentation 
for  average  noise  measurements  and  amplitude  probability  distributions  were 
included  in  the  equipment  at  each  station.  Loop  antennas  and  other  associated 
equipment  was  available  at  each  site  for  calibration  checks  and  operational 
checks  of  the  equipment.  Esterline  Angus  Recorders  were  used  to  record 
vertical  electric  field  intensities,  atmospheric  noise  levels  and  phase  variations. 
Amplitude  probability  distributions  wete  recorded  manually.  The  noise 
measuring  instrumentation  was  semi-automatic  providing  for  five  minute 
integrations  at  twelve  different  frequencies  which  recycled  automatically 


every  hour.  An  electronically  advanced  camera  ayatem  with  aaaociated 
amplifier  a  and  croeaed  loopa  waa  available  at  both  aitea  for  making  direction 
of  arrival  meaaurementa  on  the  atmoapheric  noiae. 

3.  FIELD  OPERATIONS 

The  operation  waa  beaet  with  difficultiea  and  delaya  from  the  outaet. 

The  normal  problema  of  complex  operatlona  in  the  field  were  multiplied 
by  the  tight  time  schedule,  power  aource  problema  in  both  Japan  and  Guam 
and  ahipping  damage.  The  late  arrival  of  the  equipment  at  Boulder  resulted 
in  inadequate  training  of  peraonnel  and  inadequate  laboratory  testing  of  the 
equipment  prior  to  shipment.  Some  of  the  special  equipment  such  as  the 
semi-automatic  noiae  measuring  gear  and  the  direction  of  arrival  equipment 
had  to  be  shipped  unfinished  for  completion  in  the  field.  Some  of  the  equipment 
was  temporarily  delayed  in  ahipping  overseas  and  much  of  it  arrived  in  a 
somewhat  damaged  condition. 

The  final  site  locations  selected  at  both  Guam  and  Japan  were  little 
used  airfields.  Site  location  and  quarters  in  Japan  were  excellent  and  once 
the  50  cycle  power  problem  was  compensated  for,  this  field  site  was  very 
adequate  and  provided  good  support.  The  field  site  at  Guam  required 
operation  from  a  tent  which  presented  some  operational  problems,  but  the 
most  serious  difficulty  was  involved  with  the  portable  10  kilowatt  generator 
provided  for  a  power  source.  The  equipment  proved  to  be  unreliable  and 
on  several  occasions  "ran  away"  causing  difficulties  with  the  equipment  due 
to  excessive  voltage.  Lack  of  protection  for  the  equipment  and  lack  of  a 
reliable  power  source  resulted  in  considerable  loss  in  data  since  it  was 
impossible  to  man  the  sites  twenty-four  hours  a  day  for  the  entire  measurement 
period.  In  other  words,  even  though  the  equipment  was  capable  of  recording 
data  in  an  unmanned  operation  this  could  not  be  done  because  of  possible 
damage  from  a  "run  away"  generator. 


In  part,  because  of  the  field  conditions  and  in  part  because  of 
inadequate  preparations  prior  to  going  to  the  field,  the  data  obtained  from 
the  operation  at  Guam  consisted  of  continuous  phase  data  from  both  transmitters 
and  vertical  electric  field  intensity  data  from  one  transmitter  at  a  time.  The 
vertical  electric  field  intensity  data  was  obtained  with  a  302-A  wave  analyser 
since  adequate  calibration  of  the  Tracors  at  this  site  for  amplitude  measurements 
was  not  accomplished.  The  automatic  noise  equipment  was  in  operation  for 
only  the  last  few  weeks  and  this  data  is  in  question  and  is  not  considered  reliable. 
Because  of  the  bandwidth  of  the  302-A  wave  analyser,  noise  measurements 
during  the  open  key  operation  for  the  transmitters  were  generally  too  low  to 
obtain  an  accurate  value.  The  field  intensity  data  obtained  at  Guam  appears 
to  be  very  good,  however,  and  the  phase  data  will  probably  prove  of  value  if 
completely  analyzed  at  some  later  date. 

At  the  Japan  site  some  vertical  electric  field  intensity  data  was  obtained 
at  two  frequencies  simultaneously  and  single  frequency  measurements  were 
made  for  the  greater  percentage  of  time.  Phase  measurements  at  Japan  also 
were  obtained  for  most  of  the  measurement  period.  The  automatic  noise 
equipment  at  Japan  was  put  into  operation  during  June  and  an  adequate  supply 
of  noise  data  was  obtained  at  this  site.  Amplitude  probability  distributions 
of  the  atmospheric  noise  were  also  obtained  at  this  site.  Noise  measurements 
at  the  exact  frequencies  of  transmission  were  also  impossible  at  Japan 
because  of  the  bandwidth  problems. 

Measurements  of  the  inverse  distance  field  of  both  transmitters  were 
made  in  Hawaii  as  a  means  of  determining  close  relationships  between  trans¬ 
mitter  current  and  radiated  power.  This  information  was  used  to  reduce  and 
analyse  the  data  obtained  and  is  presented  in  Appendix  B. 
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Both  noise  and  vertical  electric  field  intensity  data  were  obtained  at 
Boulder  on  a  continuous  basis  with  very  little  difficulty.  The  operation  at 
Boulder  was,  of  course,  under  laboratory  conditions  with  a  reduced  number 
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of  measurements.  No  direction  of  arrival  data  was  obtained  which  was 
considered  to  be  of  sufficient  quality  for  presentation  in  this  report. 

4.  FINAL  RESULTS 

In  spite  of  the  many  difficulties  encountered  the  program  produced 
a  large  amount  of  high  quality  data,  the  analysis  of  which,  has  proved  to  be 
most  useful  and  of  value  to  an  understanding  of  VLF  propagation  and  VLF-ELF 
atmospheric  noise.  The  data  will  be  presented  in  several  different  forms  in 
this  section  as  a  means  of  gleaning  the  maximum  amount  of  information 
from  it. 

4. 1  Continuous  Signal  Plots 

Data  from  Guam  and  Japan  have  been  plotted  together  to  provide  a 
detailed  study  of  the  difference  in  the  propagation  characteristics  along  these 
two  paths  of  almost  identical  length.  Both  paths  are  all  sea  paths,  the  only 
difference  being  in  their  direction  of  propagation.  Figure  2  presents  the 
data  obtained  at  16.  6  kHz  during  the  two,  four  day  periods  this  frequency  was 
transmitted.  Such  a  small  amount  of  data  was  obtained  that  very  little  can 
be  said  except  that  the  results  from  the  two  paths  appear  to  be  similar. 

Figure  3  presents  the  large  volume  of  data  obtained  at  19.  8  kHz. 
Although  the  mean  field  intensities  at  the  two  locations  are  very  close  to 
the  same,  a  close  examination  of  the  data  reveals  significant  differences.  It 
should  be  remembered  that  the  time  for  all  of  this  data  is  Universal  Time 
and  the  all  daylight  conditions  for  the  propagation  path  occurs  from 
approximately  2000  UT  to  0600  UT.  Although  there  is  some  fluctuation, 
there  is  generally  close  correlation  between  the  field  intensity  at  the  two 
locations  during  all  daylight  path  conditions.  It  is  noted,  however,  that 
nighttime  values  at  Japan  were  generally  2  to  5  db  higher  than  the  daytime 
values,  whereas  at  Guam  the  nighttime  field  intensities  at  this  frequency  were 
from  1  to  10  db  lower  than  the  daytime  values.  Examination  of  Figure  D-l 
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in  Appendix  D  indicates  rather  stable  conditions  should  exist  at  a  distance  p 

of  6  megameters  during  daylight  hours,  but  that  considerable  instability 

j 

might  be  expected  during  the  night.  This  is  borne  out  by  the  results  shown  !  j 

i  •  i 

in  Figure  3.  In  addition,  since  the  difference  in  distance  of  the  two  paths  '  \ 

is  less  than  1 00  kilometers  we  must  assume  that  the  difference  in  nighttime 
field  intensities  at  the  two  sites  must  result  from  different  phase  velocities 
along  the  two  paths.  Examination  of  Figures  4,  5  and  6  and  Figures  D-2, 

D-3  and  0-4  tend  to  further  support  this  analysis.  From  Figure  D-4  we  note 

i 

that  we  might  expect  greater  variations  between  Japan  and  Guam  during 

daytime  hours  than  at  night  for  26.1  kHs.  Figure  6  tends  to  bear  this  out 

since  we  observe  that  nighttime  values  at  Japan  and  Guam,  despite  their  ']  | 

scatter,  are  approximately  the  same  whereas  the  daytime  fields  at  Guam  i 

| 

are  6  to  8  db  higher  than  those  at  Japan.  Intermediate  conditions  are  seen 

to  exist  at  22.  3  and  24.  0  kHz.  * 

r,  t 

At  all  frequencies  it  is  noted  that  considerable  stability  exists  during  i 

daylight  hours.  i 

'  l  ; 
1  *  i 

Simultaneous  measurements  were  made  at  Japan  at  two  frequencies 
during  the  1 6  days  shown  in  Figure  7.  Although  normalised  absolute  field 
intensities  may  show  as  much  as  10  db  difference  between  the  two  frequencies, 
it  is  interesting  to  note  that  the  variation  in  field  intensity  at  the  two  frequencies 
is  almost  identical  with  the  exception  of  nighttime  hours  when  more  scatter 
is  seen  at  the  higher  of  the  two  frequencies.  With  the  two,  eight  day  periods 
of  data  arranged  one  above  the  other,  and  since  the  same  pairs  of  frequencies 
were  used  for  the  eight  day  periods,  it  is  possible  to  compare  results  obtained 
two  weeks  apart.  The  repeatability  of  the  results  seen  on  this  Figure  indicates 
that  with  greater  understanding  of  the  earth  ionospheric  waveguide  that 
calculations  of  expected  VLF  field  intensities  may  be  possible  even  during 
the  nighttime  and  transitional  periods.  It  appears  that  there  is  not  so  much 
a  scattering  of  data  as  there  is  a  lack  of  understanding  of  the  variations  seen. 
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4.  2  Diurnal  Signal  Plotg 


Figure*  8  through  12  are  the  plot*  of  received  signal  field  intensity 
at  Japan  aa  a  function  of  time  of  day.  It  ia  important  to  note  that  theae  value* 
are  plotted  againat  Unr-ersal  Time  and  that  local  time  ia  equal  to  UT  +  9. 

The  median  for  each  group  of  hour  and  half  hour  value*  waa  computed  and  a 
line  drawn  between  theae  median  point*.  The  xig  xag  effect  which  can  be  noted 
particularly  for  the  reaults  ai  22.  3  kHs  and  26.1  kH*  ia  the  reault  of  a  change 
in  propagation  conditiona  during  the  two  month  period  which  ahowa  up  here 
because  half  hour  values  were  obtained  only  during  the  first  month.  Similar 
results  for  Guam  are  shown  in  Figures  13  through  17  and  the  results  for  Boulder 
are  shown  in  Figures  18  through  22.  It  is  interesting  to  note  that  both  Japan 
and  Guam  experienced  changes  in  propagation  characteristics  from  the  first 
month  to  the  second  month  whereas  the  Hawaii  to  Boulder  path  doea  not  exhibit 
auch  a  change.  There  are  significant  differences  also  which  can  be  noted  between 
the  changes  observed  at  Japan  and  Gua m ,  such  as  at  22.  3  kHs  the  half  hour 
medians  plot  consistently  higher  than  the  hourly  medians  during  the  daytime 
at  both  Japan  and  Guam  whereas  during  the  night  at  Japan  the  half  hour  medians 
are  still  higher  than  the  hourly  medians,  but  at  Guam  they  are  lower.  This 
indicates  that  during  the  first  month  of  the  measurements,  propagation  was 
better  at  all  hours  at  Japan  than  during  the  second  month,  whereas  at  Guam 
propagation  during  the  first  month  was  better  during  the  day  and  worse  during 
the  night.  A  considerable  refinement  in  our  under  standing.,  of  ionospheric 
and  latitude  changes  as  a  function  of  season  will  be  required  before  our  prediction 
capabilities  can  account  for  such  variations.  Sunrise  and  sunset  times  are  at 
the  surface  and  are  for  June  15. 

4.  3  Probability  Distributions  of  Signal 

Figures  23  through  37  give  the  plots  of  the  probability  distribution  of 
the  measured  signal  intensity  at  all  three  locations,  Japan,  Guam  and  Boulder. 


Only  hour  data  was  used  to  prepare  the  probability  distribution  plots  for 
Japan  and  Guam  since  inclusion  of  half  hour  data  would  tend  to  bias  the  results 
in  favor  of  the  conditions  existing  during  the  first  month  of  the  measurement 
period.  The  generally  inverted  S  shape  of  these  curves  is  expected  because 
of  the  relative  stability  of  the  signal  during  a  large  portion  of  the  day  and  the 
considerable  instability  experienced  during  the  transition  periods  from  day 
to  night  and  night  to  day.  The  variation  of  the  signal  under  such  conditions 
is  not  random  and  we  would  not  expect  the  straight  line  plot  of  a  normal 
distribution. 

4.  4  Atmospheric  Noise  Measurements  in  Japan 

Figures  38  through  43  present  the  results  of  atmospheric  noise  measure¬ 
ments  obtained  in  Japan.  A  measurement  was  taken  approximately  once  an 
hour  at  each  frequency  for  a  period  of  31  days  from  June  14,  1965  to  July  14,  1965. 
Each  measurement  consisted  of  a  five  minute  integration  of  the  noise  in  a 
7  Hs  bandwidth.  Each  of  these  five  minute  integrations  have  been  plotted  on 
the  curves  shown.  The  data  has  been  presented  in  four  hour  time  blocks 
rather  than  for  hourly  values  for  comparison  with  CCIR  predictions  and  other 
measurements  presented  in  a  similar  manner.  The  data  has  been  normalised 
to  a  1  Hs  bandwidth  for  this  presentation. 

Continuous  noise  data  was  not  obtained  at  the  Guam  site.  The  data  which 
was  obtained  has  been  used  only  for  the  preparation  of  signal  to  noise  probability 
distributions  which  will  be  presented  in  Section  5. 

The  atmospheric  noise  data  obtained  in  Japan  at  10,  20  and  30  kHs  was 
also  plotted  versus  time  of  day  to  determine  diurnal  variations.  These  plots 
are  given  in  Figures  44  through  46.  The  reduction  in  local  thunderstorm  activity 
and  activity  in  the  Indonesian  area  is  noted  on  these  curves  during  the  hours 
2000  to  0400  UT.  A  slight  dip  can  also  be  noted  at  about  1000  UT  corresponding 
to  the  minimum  worldwide  thunderstorm  activity.  The  data  was  normalised 
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to  a  100  Hs  bandwidth  for  this  presentation. 

Atmospheric  noise  measurements  were  also  obtained  at  Boulder  and 
die  diurnal  plots  of  the  Boulder  noise  data  at  each  of  the  transmitter 
frequencies  are  shown  in  Figures  47  through  51.  It  should  be  noted  that  local 
time  in  Boulder  is  given  by  UT  -  7  hours.  A  decrease  in  local  activity  is 
noted  for  the  hours  0800  to  1 700  UT.  The  local  activity  tends  to  mask  any 
evidence  of  the  decrease  in  worldwide  activity  at  1000  UT.  The  bandwidth  is 
100  Hs. 

The  probability  distributions  of  the  rms  atmospheric  noise  field  intensities 
for  Japan  and  Boulder  are  shown  in  Figures  52  through,  67.  It  will  be  noted 
that  this  data  more  generally  approaches  a  log  normal  distribution  than  did  the 
signal  data.'  There  is  still,  however,  some  tendency  for  an  inverted  S  shape 
indicating  that  the  data  is  not  entirely  random.  This  may  be  due  tmthe  relatively 
small  number  of  data  points  available.  All  of  the  probability  plots  of  noise 
have  been  computed  from  all  data  obtained  during  the  entire  two  month  period. 

The  Japan  data  is  normalised  to  a  1  Hs  bandwidth,  Boulder  to  a  100  Hs  band¬ 
width. 

It  will  be  noticed  when  examining  the  noise  data  for  Boulder  that  there 
are  abnormal  differences  in  mean  or  median  field  intensities  for  the  relatively 
slight  difference  in  the  frequencies  measured.  This  is  logically  explained  by 
noting  that  noise  measurements  at  Boulder  were  not  made  continuously  at 
each  frequency,  but  rather  were  made  in  four  day  groups  corresponding  to 
the  transmitted  frequencies  during  which  considerable  change  in  local 
thunderstorm  activity  could  occur. 

4.5  Amplitude  Probability  Distributions 

Amplitude  Probability  Distributions  were  obtained  only  in  Japan  during 
this  program.  In  addition  to  providing  information  regarding  the  cumulative 
characteristics  of  the  atmospheric  noise  at  the  frequencies  studied,  these  APD’s 


(Amplitude  Probability  Distribution*}  provide  the  only  measurement  of  the 
average  to  rmi  conversion  factor  necessary  to  convert  the  measurements 
of  average  electric  field  intensity  to  rms  electric  field  intensity.  The  interim 
reports,  prepared  during  the  progress  of  the  measurements,  used  CCIR 
predictions  of  the  average  to  rms  conversion  factor,  since,  of  course,  the 
APD's  were  not  available  at  the  time.  Over  100  APD's  were  taken  during 
the  two  month  measurement  period  and  were  reported  in  an  interim  report. 

Two  samples  at  each  of  the  frequencies  used  are  included  in  this  final  report 
as  examples  of  low  and  wide  dynamic  range  conditions  experienced  during 
die  measurement  period.  Table  1  shows  the  minimum  and  maximum  dynamic 
range  (90  percent  to  0.  0001  percent  exceedance  level)  observed  during  the 
two  month  period  in  Japan.  The  sample  APD's  shown  in  Figures  68  through 
85  are  not  necessarily  those  that  gave  the  minimum  and  maximum  dynamic 
ranges  but  are  representative  of  the  results  obtained.  A  plot  of  rms  to  average 
conversion  factors  obtained  from  the  APD's  is  given  in  Figure  86  along  with 
values  obtained  at  other  locations  at  other  times.  The  heavy  curve  shown 
indicates  the  values  used  for  converting  the  average  noise  measurements  to 
rms  values.  Since  Malta,  Summer  and  Japan,  Summer,  conditions  are  similar, 
and  since  the  data  obtained  at  Japan  was  limited  in  quantity,  the  results  from 
both  locations  were  considered  in  constructing  the  conversion  curve. 

5.  ANALYSIS  OF  RESULTS 

In  addition  to  the  reduction  of  the  various  collections  of  signal  and  noise 
data  and  their  presentation  in  various  forms  individually,  signal -to -noise -ratios 
have  been  computed  from  the  data  and  are  presented  in  this  section.  Since  the 
variations  of  signal  and  atmospheric  noise  field  intensities  are  both  functions 
of  time  of  day  and  other  propagation  variables,  they  cannot  be  treated  as  being 
non -correlated.  It  is  important,  therefore,  when  computing  the  diurnal  pattern 
and  cumulative  distribution  of  signal-to-noise-ratios,  that  the  signal  and  noise 
field  intensities  used,  be  only  those  measured  at  or  near  the  same  real  time. 
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Since  the  Japan  noiee  data  wai  obtained  with  instrumentation  completely 
geparate  from  that  uaed  for  determining  gignal  field  intensities,  the  analysis 
of  the  data  in  terms  of  signal-to-noise-ratios  required  the  selection  of  signal 
data  and  noise  data  which  had  been  measured  within  30  minutes  of  each  other. 
Although  this  resulted  in  a  reduction  in  die  number  of  data  points  for  signal* 
to-noise- ratios  compared  to  the  total  number  of  data  for  signal  or  noise 
separately,  it  was  considered  necessary  to  analyse  die  data  in  this  manner 
to  obtain  valid  results. 

As  was  mentioned  earlier,  much  of  the  noise  data  at  Guam  did  not 
yield  absolute  field  intensities  because  of  the  dynamic  range  limitations  of 
the  field  intensity  equipment  and/or  equipment  noise  levels.  Therefore, 
diurnal  plots  of  signal-to-noise-ratios  for  Guam  are  not  available.  The 
chart  recordings  obtained,  however,  do  provide  measurements  of  the  low 
signal *10*00180 -ratios  and  also  permit  the  assignment  of  greater  than  values 
to  signaI*to*noise*ratios.  It  was  possible,  therefore,  to  manually  tabulate 
the  low  signal -to- noise -ratios  and  greater  than  values.  The  collection  of 
these  data  permitted  the  calculation  and  plotting  of  the  cumulative  distribution 
of  the  low  signal-to*noise-ratios  or  high  probability  values.  All  of  the  results 
presented  in  the  following  sub  sections  have  been  normalized  to  a  bandwidth 
of  1 00  Hs  and  a  transmitter  radiated  power  of  i  kilowatt. 

Appendix  C  describes  methods  used  to  compute  the  signal-to*noise-ratios 
described  herein. 

5.  i  Diurnal  Sgnal-to-Noise  Plots 

Figures  87  through  91  are  diurnal  plots  of  signal-to-noise-ratios 
for  Japan  and  Figures  92  through  96  are  diurnal  plots  of  signal-to-noise-ratios 
as  measured  at  Boulder.  Sunrise  and  sunset  times  were  computed  at  the 
earth's  surface  on  June  15. 

An  examination  of  Figures  87  through  91  shows  very  similar  diurnal 
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variation*  of  the  signal-to-noise-ratios  with  the  exception  of  16.  6  kHa  which 
indicate*  a  generally  poor  signal-to-noise-ratio  from  aonaet  to  aunriae.  The 
reault*  for  19.  8.  22.  3  and  24.  0  kHa  are  eapecially  aimilar  and  in  thia  normalised 
form  indicate  that  for  a  given  power  radiated  equivalent  communication*  might 
be  expected  at  these  frequencies.  Reference  to  Appendix  B,  however,  which 
shows  a  general  increase  in  radiated  power  with  an  increase  in  frequency  would 
indicate  that  better  results  might  be  obtained  at  the  two  higher  frequencies.  At 
26.1  kHs,  generally  lower  normalised  signal -to -noise -ratio*  are  indicated. 

When  examining  Figures  92  through  96  for  Boulder  it  should  be  noted 
that  local  time  is  now  UT  -  7  not  UT  +  9  as  was  the  case  for  the  Japan  data. 

The  data  for  16.  6  and  19.  8  kHs  show  a  reduced  diurnal  variation  with  a  broad 
sunset  effect  from  0200  to  0500  UT  and  a  short  sunrise  dip  at  1500  UT.  The 
higher  frequencies  show  an  increasingly  greater  reduction  in  signal  during  the 
sunrise  and  sunset  periods  with  markedly  higher  signal -to-noise -ratios  during 
nighttime  hours.  The  sunset  effect  is  consistently  spread  over  a  number  of 
hours  from  0200  to  0500  with  the  minimum  signal-to-noise-ratio  changing  with 
frequency  from  0330  at  19.  8  to  0500  at  22.  3  kHs.  The  sunrise  effect  is  the 
same  at  all  frequencies  with  a  sharp  dip  in  signal-to-noise-ratio  occurring 
consistently  at  1500  UT.  The  saw  tooth  effect  or  peaks  or  valleys  occurring 
on  the  half  hour  and  hour  for  the  plot  of  median  signal-to-noise-ratios  is  in 
this  case  an  indication  of  a  change  in  atmospheric  noise  conditions  from  the 
first  half  of  this  two  month  period  relative  to  the  second  half  since  such  a 
sig  sag  effect  was  not  noted  for  signal  data  at  Boulder.  Increased  local 
thunderstorm  activity  for  the  second  30  day  period  resulted  in  a  lower  signal- 
to-noise-ratio  as  expected. 

5.  2  Probability  Distributions  of  Signal-to-Noise -Ratio 

Figures  97  through  101  are  plots  of  the  cumulative  distribution  of 
measured  signal-to-noise-ratios  for  Japan,  Figures  102  through  108  are  similar 
figures  for  Guam  and  Figures  109  through  113  are  for  Boulder.  It  is  important 
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to  remember  that  the  Japan  data  wai  obtained  by  interpolation  of  noiee  data 
measured  at  1 0,  20  and  30  kHa  to  the  frequencies  of  interest  and  by  selecting 
signal  and  noise  data  measured  during  the  same  half  hour  period.  Figures  102 
through  108  for  Guam  were  manually  computed  as  previously  described  and  it 
should  be  noted  that  data  for  Lualualei  and  Haiku  have  been  presented  separately 
to  provide  a  comparison  of  the  results  with  this  technique.  It  is  encouraging 
to  note  that  Lualualei  and  Haiku  are  in  close  agreement.  The  data  for  Boulder 
includes  all  data  collected  during  the  two  month  period  and  it  is  noted  that 
with  this  very  much  greater  number  of  data  points,  approximately  10  times 
the  data  available  from  Japan,  that  the  plots  are  very  nearly  log  normal. 

These  plots  permit  the  determination  of  the  percentage  of  time  that 
a  communications  system  requiring  a  given  signal-to-noise-ratio  will  be 
effective  in  these  areas. 

Complete  details  on  the  methods  of  computation  for  these  cumulative 
distributions  is  given  in  Appendix  C. 

5.  3  Means,  Medians  and  Standard  Deviations 

Tables  2  through  12  give  the  computed  values  for  the  mean,  median  and 
standard  deviation  of  signal,  noise  and  signal-to-noise-ratio  for  Japan,  Guam 
and  Boulder.  Four  hour  time  block  presentations  are  included  along  with  the 
values  for  all  twenty-four  hours  combined.  There  were  no  computations  for 
noise  and  signal-to-noise-ratio  for  Guam  since  the  noise  data  collected  at 
that  location  was  incomplete.  For  Japan  and  Guam  the  haLi  hour  data  was 
omitted  from  computations  for  the  signal  values  since  its  inclusion  would  have 
weighted  the  results  in  favor  of  the  first  half  of  the  two  month  period.  The 
Japan  noise  and  signal -to -noise  computations  were,  of  course,  the  selected 
data  as  previously  described  which  covered  a  31  day  period.  All  data  collected 
at  Boulder  was  used  in  the  Boulder  computations. 

Tables  13  and  14  give  the  computations  of  standard  deviation,  mean 
and  median  values  for  the  atmospheric  noise  data  obtained  in  Japan  at  frequencies 
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from  20  Ha  to  30  kHs.  In  contrast  to  the  noise  data  used  to  obtain  signal- 
to- noise -ratios,  this  data  has  been  normalised  to  a  i  Hs  bandwidth  instead 
of  a  100  Hs  bandwidth. 

The  equations  and  methods  used  in  computing  the  standard  deviation, 
mean  and  medians  as  given  in  these  tables  are  described  in  detail  in 
Appendix  C. 

6.  RESULTS  OF  PHASE  MEASUREMENTS 

The  analysis  of  the  phase  data  to  date  has  been  very  time  consuming 
and  disappointing.  Phase  jumps  and  excessive  phase  drift  of  the  transmitter 
frequency  standard,  the  receiver  frequency  standard  or  both  has  rendered 
much  of  the  data  extremely  difficult,  if  not  impossible  to  analyse.  In  addition, 
we  find  numerous  occasions  when  the  Tracor  phase  receivers  lost  tract  or 
were  shut  down  for  one  reason  or  another,  including  primary  power  failures. 
Given  sufficient  time  to  compare  results  between  measurements  made  at 
the  same  time  with  different  pieces  of  equipment  and  to  correlate  phase  dis¬ 
continuities  with  transmitter  records,  a  large  quantity  of  good  and  useful 
phase  data  could  be  obtained  from  the  raw  data  in  hand.  Since  it  was  considered 
that  the  phase  data  was  of  lesser  importance  than  the  amplitude  and  signal-to- 
noise-ratios,  the  analysis  of  the  phase  data  was  left  until  last.  Time  and  fiscal 
limitations  and  the  difficulties  described  above  have  limited  the  amount  of 
data  available  for  presentation  in  this  report. 

Figures  114  through  118  present  the  results  of  the  data  analyzed  at 
this  time.  Figures  114  and  115  show  the  diurnal  phase  variations  as  measured 
at  Japan  and  Guam  during  the  four  day  period  from  26  June  to  29  June  1965. 

The  phase  variations  at  Japan  are  reversed  from  those  at  Guam  and  are 
reversed  from  those  predicted  by  theory.  Although  this  could  have  been  caused 
by  improper  connections  of  the  equipment  there  is  no  record  to  indicate  that 
such  was  the  case.  We  are  not  discounting  the  fact  that  this  could  have 
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happened,  however.  The  discontinuity  indicated  on  Figure  115  indicates  a 
phase  jump  which  could  not  be  accurately  determined.  The  phase  plot  from 
0800  to  2400  on  the  28th  of  June  was  arbitrarily  positioned  such  that  die  relative 
phase  observed  at  2400  hours  on  the  27th  of  June  was  identical  to  that  at  2400 
hours  on  the  28th  of  June. 

Figure  116  shows  the  diurnal  variation  of  phase  occurring  on  the 
17th  of  June  at  a  frequency  of  22.  3  kHs  in  Japan.  Assuming  die  plotting  of 
the  phase  variations  following  the  discontinuity  has  been  accurately  positioned, 
it  appears  that  phase  lags  occur  at  sunrise  and  sunset  but  the  phase  during 
the  nighttime  and  daytime  periods  is  approximately  the  same. 

Figures  117  and  118  present  the  diurnal  variation  of  phase  at  26.1  kHz 
for  the  period  22  June  to  25  June.  A  much  greater  diurnal  variation  is  noted 
at  Japan  than  at  Guam.  It  is  particularly  interesting  to  note  that  this  four 
day  period  immediately  preceding  the  four  day  period  presented  in  Figure  114  and 
that  a  normal  diurnal  variation  of  phase  was  obtained  at  26. 1  kHz  whereas 
a  reversed  diurnal  variation  was  obtained  at  19.  8  kHs.  There  is  no  reason 
to  believe  that  any  changes  in  the  equipment  connections  were  made  betw  an 
these  two  four  day  records.  Further  comparisons  of  Figures  114,  116  and 
117  would  lead  one  to  believe  that  at  19.  8  kHz  in  Japan  a  reversed  phase 
variation  was  obtained,  at  22.  3  kHz  no  significant  nighttime  to  daytime 
variation  in  phase  was  observed  and  that  at  26. 1  kHz  a  significant  and  normal 
phase  variation  was  observed. 

As  mentioned  above  the  examination  of  considerable  additional  data 
should  be  done  before  any  conclusions  are  reached  concerning  the  diurnal 
phase  variations  at  these  locations.  , 

7.  CONCLUSIONS 

The  primary  conclusion  which  this  data  leads  one  to  make  is  simply 
that  the  state-of-the-art  today  is  not  adequate  for  accurately  predicting  diurnal 
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variation*  of  VLF  field  intensities.  In  fact,  it  is  barely  accurate  for  the 
prediction  of  monthly  means  and  medians.  This  is  especially  true  at  frequencies 
above  20  kHs. 

This  measurement  program  has  resulted  in  a  significantly  greater 
understanding  of  field  intensities  and  signal-to-noise -ratios  in  the  western 
Pacific  region.  For  communications  between  NPM  and  this  western  Pacific 
area,  either  of  the  frequencies  22.  3  or  24.  0  kHz  should  result  in  significantly 
improved  communications  during  this  season  than  16.6,  19.8  or  26.1  kHs. 

It  should  not  be  assumed  that  this  will  be  true  for  other  seasons. 

The  opportunity  to  simultaneously  measure  field  intensities  over 
identical  or  similar  paths  at  identical  frequencies  and  to  make  measurements 
at  different  frequencies  over  the  same  path  simultaneously  has  proven  to  be 
of  significant  value.  The  repeatability  of  diurnal  variations  over  long  periods 
of  time  and  some  significant  differences  noted  over  very  similar  paths  as  a 
function  of  frequency  provide  an  unparalleled  opportunity  for  improving  our 
understanding  of  VLF  propagation. 


It  is  interesting  to  note  that  a  change  in  procedure  (eliminating  half 
hour  measurements)  from  the  urst  30  day  period  to  the  second  30  day  period 
accentuated  changes  in  propagation  characteristics  and  atmospheric  noise 
levels.  These  differences  might  never  have  been  noted  without  this  fortuitous 
change  in  operation.  This  should  be  kept  in  mind  for  future  programs. 

Finally  it  can  be  stated  without  reservation  that,  although  very 
significant  results  have  been  forthcoming  from  this  program,  that  they  could 
have  been  obtained  with  much  less  difficulty  and  in  the  long  run,  in  less 
time  had  the  operational  schedule  been  planned  more  realistically  or  the 
number  of  measurement  programs  reduced  significantly.  This  too  should 
provide  guidance  for  future  programs. 

It  is  recommended  that  further  analysis  of  this  data  be  made,  in 
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particular,  the  phase  data  available  and  that  comparisons  of  the  data  be 
made  with  more  rigorous  computations  of  field  intensity  including  the  effects 
of  mode  conversion  at  the  sunrise-sunset  interfaces.  It  is  obvious  that  some 
variation  in  the  characteristics  of  the  ionosphere  for  the  paths  must  also  be 
considered. 
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TABLE  i 


Dynamic  Range  (90f.  to  0.  0001  %  Exceedance  Levels)  of 
Atmospheric  Noise  Observed  in  Japan. 


Freq.  Hx 

Minimum 

Maximum 

140 

47 

61 

no 

39 

61 

•520 

43 

68 

1  K 

43 

68 

3  K 

36 

62 

5  K 

36 

72 

10  K 

37 

73 

20  K 

41 

67 

30  K 

41 

60 
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Standard  Deviation  in  db 


Standard  Deviation,  Mean  and  Median  of  Selected  Noiae 
From  Japan  For  Four  Hour  Time  Blocks. 
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Mean/ Median  in  db. 
Standard  Deviation  in  db 


TABLE  7 

Standard  Deviation,  Mean  and  Median  of  Hourly  Signal 
From  Guam  For  Four  Hour  Time  Blocka. 


Standard  Deviation,  Mean  and  Median  of  all  Signal  Data 
From  Boulder  For  Four  Hour  Time  Rinrk. 
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Figure  9  Diurnal  Variation  of  Signol 
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Figure  13  Diurnal  Variation  of  Signal 
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t  17  Diurnal  Variation  of  Signal 
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Figure  18  Diumol  Variation  of  Signal 


Figure  19  Diurnal  Variation  of  Signal 


Figur*  21  Diurnal  Variation  of  Signal 
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Figure  23  Probability  Diitribution  of  Meatursd  Signal  Inteniily 


Probability  That  the  field  V.'ill  Exceed  the  Ordinate,  Percent  of  Time  99.9 

Figure  25  Prooability  Oiitribution  of  Measured  Signol  Inteniity 
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Figure  27  Probobility  Distribution  of  Measured  Signal  Intensity 


-  Summer  19 
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Probability  That  the  Field  Will  Exceed  the  Ordinate,  Percent  of  Time 
Figure  31  Probability  Distribution  of  Measured  Signal  Intensity 
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Figur#  33  Probability  OUtribution  of  M#atur#d  Signal  Intanthy 
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Figaro  35  Probability  Diitribotion  of  Maamrod  Signal  Intensity 
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Location  -  Boulder 
S«o ion  -  Summer  19 
Frequency  *  26. IK) 
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Figure  37  Probability  Distribution  of  Measured  Signal  Intensity 


Frequency  in  Hz 

Figure  4)  Atmospheric  Noise  Spectra 


-  Summar  1965 
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Figur*  44  Diurnal  Variation*  of  VLF  Aimotpharic  Noito 


Figure  45  Diurnal  Variation*  of  VLF  Atmospheric  Noise 
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Figure  46  Diurnal  Variations  of  VLF  Atmospheric  Noise 


Flgur*  47  Diurnal  Variation  of  VLf  Atmospheric  No  Iso 
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Figura  53  Cumulatlva  Distribution  of  RMS  Atmospharlc  No  is*  Intansity 


1 


8  R  $  S  9 

i/3^/Ui/Af<  1  04  SAIJOjag  qp  uj  X4|$U84U|  p|»lj 


Figure  61  Cumulative  Distribution  of  RMS  Atmospheric  Noise  Intensity 
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igure  64  Cumulative  Distribution 


Distribution  of  Meosured  RMS  Atmospheric  Noise 


Location  -  Boulder 
Season  -  Summer  1965 


Figure  66  Cumulative  Distribution  of  Measured  RMS  Atmospheric  Noise 
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Figure  67  Cumulative  Distribution  of  Measured  RMS  Atmospheric  Noise 


Vertical  Electric  Field  Intensity,  db  relative  to  1  pv/m  30-935C 

Figure  69  Amplitude  Probability  Distribution  for  Atmospheric  Noise 
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Vertical  Electric  Field  intensity,  db  relative  fa  1  pv/m 
Figure  73  Amplitude  Probability  Distribution  for  Atmospheric  Noise 
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Vertical  Electric  Field  Intensity,  db  relative  to  1  pv/m 
Figure  77  Amplitude  Probability  Distribution  for  Atmospheric  Noise 


Figure  78  Amplitude  fcobobiHty  Disfr Outran 


V^Tticol  Electric  Field  Intensity,  <ft>  relotive  to  1  u*/r» 
Figure  TV  Amplitude  Probability  Distribution  for  Atmospheric  Noil 
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Vertical  Electric  Field  Intensity,  db  relative  to  1  pv/ro 
’  82  Amplitude  Probability  Distribution  for  Atmospheric  Noi 


Vertical  Electric  Field  intensity,  db  relative  to  1  pv/m 
i’4  Amplitude  Probability  Distribution  for  Atmospheric  Noi 
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igure  88  Meosured 


[gure  89  Measured  Signal  to  Noise  Ratio 
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Figure  90  Meosured  Signal  to  Noise  Fcfio 


Figure  91  Measured  Signal  to  Noi: 
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Figure  9'J  Measured  Signal  to  Noise  Ratio 
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June  30  -  July  4 
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Figure  99  Cumulative  Distribution  of  Measured  Sign  j)  to  Noise  Rotiu 
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Figure  101  Cumulative  Distribution  of  Measured  Signal  to  Noise  Ratio 


Figure  102  Cumulative  Distribution  of  Measured  Signal  to  Noise  Ratio 


Figure  104  CuraiiloiN*  Distribution  of  Measured  5igac!  to  Noise  Ratio 


Location  -  Guam 
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Figure  109  Cumulative  Distribution  Measured  Signal  to  Noise  Ratio 


Figure  110  Cumulative  Distribution  of  Measured  Signal  to  Noise  Ratio 
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Location  “  Japan 
Seoson  -  Summer  1965 
Frequency  "  19.8  KHz 
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Locotioo  -  Japan 


igure  117  Diurnol  Variofion  of  Phase 


APPENDIX  A 

INTERIM  REPORT  ERRATA 


During  the  final  analysis  of  the  Pacific  test  data,  a  number  of  errors 
were  discovered  in  the  preliminary  analysis.  Compiled  herein  is  an  errata 
list  which  may  be  used  to  correct  the  interim  reports  (released  during  the 
summer  of  1965). 


1  boulder 

1.  All  of  the  signal  and  noise  data  must  be  raised  by  1.5  db. 

2.  June  6  through  June  9  (26.1  kHs)  -  subtract  9.  7  db. 

H  GUAM 

Due  to  the  narrow  bandwidth  of  the  receiver  used  to  measure  field 
intensities  of  the  signal  and  atmospheric  noise  field  intensities  at  Guam, 
much  of  the  noise  data  obtained  at  that  site  is  inaccurate.  It  is  suggested 
that  all  of  the  atmospheric  noise  levels  for  Guam  which  were  given  in  the 
interim  reports  be  considered  suspect.  The  data  which  was  accurate  or 
suitable  for  correction  has  been  re -tabulated  and  will  be  available  along  with 
all  other  tabulated  data  from  the  Bureau  of  Slips  and  the  Naval  Research 
Laboratory. 


A.  Data  Changes 


1. 

2. 

3. 

4. 

5. 

6. 
7. 


8. 


May  i  7  through  May  20  (24.  0  kHs)  -  subtract  20  db; 

May  22  0100  hours  (26. 1  kHs)  -  subtract  1. 5  db; 

May  24  2100  hours  (26.1  kHs)  -  subtract  10  db; 

May  29  2330  through  May  30  0300  hours  (22.  3  kHs)  -  add  10 

June  5  1630  through  1930  hours  (24.  0  kHs)  -  add  10  db; 

June  8  2300  through  2400  hours  (26.1  kHs)  -  add  10  db; 

June  10  2200  hours  (19.  8  kHs)  -  add  4.  5  db; 

June  14  2200  hours  (22.  3  kHs)  -  add  2.  0  db; 


9. 

June  15 

2200  hours  (22.  3  kHx)  -  add  3.  5  db; 

10. 

July  4 

2300  hours  (16.  6  kHx)  -  add  2.  0  db; 

11. 

July  11 

2200  hours  (19.  8  kHs)  -  subtract  6.  0  db. 

Data  to  be  Deleted 

1. 

May  18 

2300  hours  (24.  0  kHx); 

2. 

May  21 

2400  hours  (26. 1  kHx); 

3. 

May  24 

0030  through  1030  and  2130  hours  (26.1  kHx); 

4. 

May  29 

2200  and  2230  hours  (22.  3  kHx); 

5. 

May  30 

2300  hours  (22.  3  kHx); 

6. 

June  4 

2230  through  June  5  0100  hours  (24.  0  kHx); 

7. 

June  5 

0300  and  0330  hours  (24.  0  kHx); 

8. 

June  5 

2000  through  2400  hours  (24.  0  kHs); 

9. 

June  4 

1800  through  June  5  0400  hours  (19.  8  kHs); 

10. 

June  7 

0100  hours  (26.1  kHx); 

11. 

June  28 

0400  hours  (19.  8  kHx). 

m  JAPAN 

All  of  the  Japan  noise  data  except  the  Amplitude  Probability  Distributions 
as  presented  in  the  interim  reports  must  be  lowered  by  10. 1  db.  The  tables 
of  signal  data  for  Japan  in  the  interim  reports  did  not  contain  the  measurements 
made  on  the  half  hour  period.  These  tables  have  been  re -tabulated  and  will 
be  available,  if  needed  from  the  Bureau  of  Ships  and  the  Naval  Research 
Laboratory. 

A.  Data  Changes 

1.  May  17  through  May  20  (24.  0  kHx)  -  add  3.  0  db; 

2.  May  17  through  May  20  {19.  8  kHs)  -  add  4.  6  db; 

3.  May  21  through  May  24  (26. 1  kHx)  -  add  2.  4  db; 

4.  May  21  through  May  27  (19.  8  kHx)  -  subtract  9.4  db; 

5.  May  25  through  May  28  (24.  0  kHx)  -  add  3.  0  db; 
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6.  May  29  through  May  31  0600  hours  (22.  3  kHs)  -  add  3.  6  db; 

7.  Juno  1  0300  hours  (22.  3  kHs)  -  add  1 0.  0  db; 

8.  Juno  2  0700  through  June  3  1 800  hours  (1 9.  8  kHs)  -  add  3.  0  db; 

9.  June  8  0900  through  June  9  2400  hours  (1 9.  9  kHs)  -  subtract  8.  0  db; 

10.  June  13  0100  through  June  15  1200  hours  (16.6  kHs)  -  subtract  10.0  db; 

11.  June  18  0200  through  0700  hours  (24.  0  kHs)  -  subtract  10.  0  db; 

12.  June  18  0200  and  0900  hours  (24.  0  kHs)  -  subtract  20.  0  db; 

13.  June  22  1900  through  June  24  1000  hours  (26.1  kHs)  -  subtract  10.  Odb; 

14.  July  2  0000  hours  (22.  3  kHs)  -  add  10.  0  db. 

B.  Data  to  be  Deleted 

1.  May  28  0030  through  June  2  0700  hours  (19.  8  kHs). 
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APPENDIX  B 
RADIATION  RESISTANCE 
Roland  D.  Croghan 

Daring  the  period  of  time  between  June  4  and  June  10,  1965,  Reid 
strength  measurements  were  taken  on  the  transmitter  signals  from  die 
Loaloalei  and  Haiku  Naval  VLF  Transmitters  located  on  the  Island  of  Oahu, 
Hawaii.  The  purpose  of  the  measurements  was  to  calculate  the  radiation 
resistance,  as  a  function  of  frequency,  of  both  the  Lualualei  and  Haiku 
transmitters.  All  field  strength  measurements  were  taken  so  that  the  distance 

N 

from  the  transmitter,  d,  was  \/Z  <.  d  <  52  km  (this  range  is  within  the  distance 
where  the  radiation  fields  vary  inversely  with  distance  over  homogeneous 
ground  of  expected  worst-case  conductivity). 

The  actual  sites  of  measurement  are  shown  on  Figure  B-l.  The 
majority  of  the  sites  shown  were  measuring  sites  used  by  NRL  personnel 
in  previous  radiation  resistance  measurements.  DECO  Electronics  duplicated 
the  majority  of  these  measuring  locations;  in  addition,  three  new  measuring 
sites  were  added. 

During  the  measurements  Haiku  operated  on  the  frequencies  of  16.  6 
and  19.  8  kHs.  During  normal  operation  a  CW  signal  was  transmitted  at 
reduced  power.  At  scheduled  times  the  power  was  increased  for  three  minutes 
followed  by  three  minutes  of  off-time.  Lualualei  transmitted  during  this 
period  on  frequencies  of  19.  8,  24.  0  and  26.1  kHz.  Transmissions  during 
normal  operation  were  both  ICW  and  FSK.  At  scheduled  times  a  three  minute 
CW  signal  was  transmitted  followed  by  three  minutes  of  off-time. 

The  block  diagram  in  Figure  B-2  shows  the  measuring  equipment 
used  to  measure  the  radiated  component  of  the  magnetic  field.  A  loop  was 
used  with  a  calibration  box  followed  by  a  tuning  and  amplifier  section  to 
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discriminate  between  the  transmitted  signals  of  Lualualei  and  Haiku. 

Additional  tuning  vni  provided  in  the  detection  unit  along  with  two  method* 
of  detection.  Be*ide*  a  direct  reading  voltmeter,  audio  elide -back  detection 
wa*  available. 

All  measurement*  except  two  were  taken  when  the  transmitter*  were 
on  CW  transmission  (the  two  were  when  Lualualei  was  on  FSC).  After 
tuning  on  the  signal  and  maximising  the  signal  with  the  loop,  an  output  voltage 
was  read.  During  the  off  period  a  calibrated  signal  in  db  relative  to  1  pv/ m 
from  the  oscillator  was  injected  to  match  the  voltmeter  reading.  The  resultant 
values  of  field  strength  are  shown  in  Tables  1  and  11  respectively,  for 
Lualualei  and  Haiku.  The  field  strengths  have  been  normalised  to  the  indicated 
antenna  currents.  The  data  is  also  shown  in  Figures  B-3  and  B-4  where  the 
field  strength  is  plotted  as  a  function  of  distance  and  frequency.  The  spread 
in  the  data  as  shown  on  the  figures  from  an  inverse  distance  line  can  be 
contibuted  to  the  nonhomogeneity  of  the  ground,  causing  ground  impedance 
discontinuities  which  affects  the  magnitude  of  the  field.  The  best  inverse 
distance  line  fit  to  the  data  can  be  used  with  the  following  equation  to  derive 
the  radiated  power  for  the  reference  antenna  current; 

Py  J  db  rel.  1  kwj  *  J  db  rel.  i  ^v/m  j  -  109.  5  f  20  log  d  J  kmj  .  (B 
The  radiation  resistance  is  then 


1  p 

[  OhmS  j  a 

watts 

amperes2 

where  I  is  the  reference  base  antenna  current.  The  resulting  values  of 
radiation  resistance  are  shown  in  Table  HI. 


Figure  8-3  Normalized  Field  Strengths  of  the  Haiku  Transmitter 
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It  should  be  noted  that  very  few  data  points  are  available  for  24  kHs( 
so  Rf  should  not  be  weighted  as  heavy  as  the  19.  8  and  26.1  kHs  data.  The 
above  data  were  supplied  to  NRL  to  average  with  their  measurements  for 
determining  the  average  effective  height.  The  effective  height  of  the  antennas 
as  a  function  of  radiation  resistance  is: 


(B-3) 


where  h^  is  the  effective  height  in  meters  and  is  the  wavelength  in  meters. 


TAB1J5  HI 


Station 

Frequency 

kHs 

Radiation  Resistance 
ohms 

Lualualei 

19.  8 

0.  068 

24.0 

0.  082 

26.1 

0.0914 

Haiku 

16.6 

0.119 

19.  8 

0.161 

APPENDIX  C 

METHODS  OF  DATA  REDUCTION 
Donald  L.  Stone 

The  initial  data  from  each  of  the  three  receiver  station#  was  reduced 
differently  due  to  the  difference  in  the  receiving  equipment  used.  The  signal 
data  taken  at  Boulder,  using  the  URM-6,  was  contaminated  with  atmospheric 
noise  when  the  signal-to-noise-ratio  was  low.  The  Tracor  receiver  which  was 
used  in  Japan  had  a  bandwidth  narrow  enough  that  an  uncontaminated  signal 
amplitude  was  measured.  The  noise  data  from  Japan  was  obtained  with  a 
HP  302-A  receiver  in  connection  with  an  automatic  frequency  selecting 
instrument.  The  data  taken  in  Guam  with  the  HP  302-A  was  similar  to  the  URM-6 
data  and  was  reduced  in  the  same  manner. 

All  of  the  scaled  data  was  punched  onto  tabulating  cards  for  computer 
processing.  All  of  the  field  intensities,  means,  medians,  standard  deviations 
and  probabilities  were  computed  by  a  CDC  3600  digital  computer  and  plotted 
on  a  EAI  3440  Data  Plotter  with  the  exception  of  the  Guam  signal-to-noise-ratios 
and  the  phase  plots.  The  formulas  used  are  discussed  in  the  following  sections. 

SIGNAL 

A11  data  was  scaled  from  the  EA  strip  charts  in  db  relative  to  1  micro¬ 
volt.  The  contaminated  signal  data  was  first  converted  to  microvolts  and 
the  signal  computed  by  the  following  formula: 
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E  2  -  E  2 
an  n 


(C-i) 


where: 


E(  is  the  signal  in  pv. 
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E^o  is  the  contaminated  signal  in  pv,  and 

E  is  the 'noise  in  uv. 
n 

The  signal  was  then  converted  to  db  relative  to  1  pv  and  the  field  intensity 
computed  by  the  following: 


E  '  a  E  +K 

s  si 


(C-2) 


where: 

I 

E^  is  the  signal  field  intensity  in  db  relative  to  1  pv/m, 

E^  is  the  signal  in  db  relative  to  1  pv,  and 

Kj  is  a  constant  including  effective  height. 

NOISE 

The  noise  field  intensity  was  computed  by  the  following: 


E  *  E  +  K,  (C-3) 

n  n  Z 

where: 

t 

E  is  the  noise  field  intensity  in  db  realtive  toi  pv/m  for  a 
n  specified  bandwidth, 

Eq  is  the  average  noise  in  db  relative  to  i  pv,  and 

1*  a  constant  which  included  the  effective  height,  bandwidth 
factors,  and  average  to  rms  conversion  factor. 

SIGNAL-TO-NOISE  -  RA  TIO 


The  signal»to-noise-ratios  for  Boulder  and  Guam  were  computed 
by  the  following: 


t 


(C-4) 


S/N  *  E  '  -  E 

s  n 


where: 

S/N  is  the  signal -to-noise- ratio  in  db  for  the  specified  bandwidth. 

Since  the  automatic  noise  equipment  at  Japan  did  not  measure  the 
atmospheric  noise  at  the  same  frequencies  as  the  signal  transmissions, 
an  interpolated  noise  level  was  used  to  compute  the  signal-to-noise-ratios. 
For  example,  the  noise  level  for  16.  6  kHz  was  linearly  interpolated  from  the 
values  of  noise  measured  at  1 0  kHz  and  20  kHs.  To  minimize  the  error  of 
this  method,  the  only  noise  data  used  was  that  taken  within  i/2  hour  of  the 
particular  signal  measurement  and  a  separate  interpolated  noise  value  was 
used  for  each  signal  measurement. 

MEDIANS 

To  compute  the  median  of  an  array  of  data  the  array  was  first 
numerically  ordered.  (An  array  of  data  was  defined  as  any  group  of  signal 
data,  noise  data,  etc. )  The  median  value  of  the  array  was  then  defined  as 
follows,  if  N,  the  number  of  data,  was  odd 


median  =  X 


(N  +  i)/2 


however  if  N  was  even 


median  = 


XN/2  *  XN/2  +  1 
2 


(C-5) 


(C-6) 


where : 

X  is  an  element  of  the  array. 


C-3 


STANDARD  DEVIATIONS 


To  compute  the  standard  deviation  of  an  array  of  data  the  arithmetic 
mean  of  the  voltage  (not  the  log  of  the  voltage)  must  be  computed  using  the 
equation 


X 


N 


1  *  i 


(C-7) 


where: 

X  is  the  mean, 

X^  is  an  element  of  the  array, 

N  is  the  total  number  of  elements, 
i  «  1 ,  2,  3 . N. 

The  standard  deviation  was  then  computed  using 

N 

9  3  £  (X  -  X)2  (C-8) 

i  *  i 

where: 

9  is  the  standard  deviation,  and  the  other  symbols  as  defined  above 
expressed  in  db. 

CUMULATIVE  PROBABILITY  DISTRIBUTIONS 

To  compute  the  probability  that  a  given  level  was  exceeded  in  an 
array  of  data,  the  array  was  first  numerically  ordered.  The  probability  for 
a  given  level  was  then  computed  using 


100  (N  -  K.) 
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where: 

ie  the  probability,  in  percent,  for  the  ith  db  level, 

ie  the  number  of  elements  of  the  array  that  have  a  magnitude 
leee  than  the  ith  db  level,  and 

N  ie  the  total  number  of  elements  in  the  array. 


Since  the  noise  measurements  in  Guam  were  generally  inaccurate 
for  large  S/N  ratios,  it  was  impossible  to  accurately  establish  signal-to- 
noise-ratios  greater  than  20  db  in  a  7  Hs  bandwidth.  For  ratios  less  than 
20  db,  cumulative  probability  distributions  were  manually  computed  and 
plotted  using  equation  C-9. 


One  of  the  more  important  aspects  of  the  data  reduction  was  the 
elimination  of  bad  data.  Any  data  or  calibration  which  was  suspect 
for  any  reason  was  thrown  out.  Calibrations  were  made  in  the  field  with 
several  antennas  and  instruments  a  number  of  times  during  the  measurement 
period.  Some  of  the  equipment  was  also  calibrated  upon  its  return  to  Boulder. 
The  excellent  agreement  of  these  many  checks,  and  the  elimination  of  data 
when  there  was  not  agreement,  has  given  us  high  confidence  in  the  results 
from  this  program. 
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APPENDIX  D 

CALCULATION  OF  FIELD  STRENGTH  VS.  DISTANCE 


Garner  and  Rhoads  of  NRL  computed  the  variation  of  field  strength 
with  distance  from  0.  5  to  10.  0  megameters,  taking  into  account  the  1st,  2nd 
and  3rd  order  waveguide  modes  for  all  distances.  This  was  done  for  ionosphere 
heights  of  70  and  90  km.  Figures  D-l  to  D-4  present  the  results  of  their 
computations  for  the  frequencies  19.  8,  22.  3,  24.  0  and  26. 1  kHs.  These 
figures  are  included  for  comparison  with  the  measured  fields  and  are  discussed 
in  Section  4. 1 . 


Figure  D-2  Computed  Field  Strength 
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